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Zh ' In this paper, wc calculate the stable Wigner-Seitz (W-S) cells in the inner crust of 

i ' neutron stars and we discuss the nuclear shell effects. A distinction is done between the 

shell effects due to the bound states and those induced by the unbound states, which 
are shown to be spurious. We then estimate the effects of the spurious shells on the total 
' energy and decompose it into a smooth and a residual part. We propose a correction to 

\J^i ' the Hartree-Fock binding energy in Wigner-Seitz cell (HF-WS). 
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1. Introduction 



[■ — ' From the very first nuclear models based on extrapolations of a mass formula, 1 

several models have been developed to provide the equation of state in the inner 
crust, like the BBP 2 and BPS 3 semi-classical models or fully quantal self-consistent 
mean-field models. 4 ~ 7 ' 9 ' 10 The latter models are based on the Wigner-Seitz (W-S) 
approximation and the continuum states are discretized in a box (spherical or not). 
In semi-classical models, a clear limitation is coming from the absence of nuclear 
shell effects which is known to play an important role in determining the compo- 
sition of the ground state prior to the neutron drip, 4 in the inner crust 6 as well 
as in the transition zone between the crust and the core. 11 Since the seminal work 
of Negele & Vautherin, 4 many efforts have been invested in the developpement 
of self-consistent mean- field modelization of the inner crust, by introducing defor- 
mation of the nuclear cluster in the high density part of the crust, 5,10 or pairing 
correlations, 7-9 or more recently by improving of the lattice description within the 
band theory. 12 Band theory takes into account the proper symmetries of the sys- 
tem, but the equations are numerically very complicated to solve and have not yet 
been solved in a self-consistent framework. Instead, it has been coupled to a self- 
consistent mean- field model at the W-S approximation. 13 It has then been shown 
that the W-S approximation is justified if the temperature is larger than about 



February 2, 2008 6:48 



WSPC - Proceedings Trim Size: 9.75in x 6.5in margueron'exo 



2 



Table 1. Comparison of the proton number Z and the radius of the W-S cell 
Rws obtained in different calculations. See the text for more explanations. 



zone 



kp Z Rws ( m f m ) 

fm" 1 N-V 4 HF-BCS a HF HFB N-V 4 HF-BCS a HF HFB 



2 
3 
4 
5 
6 
7 
8 
9 

10 
11 



1.12 
0.84 
0.64 
0.55 
0.48 
0.36 
0.30 
0.26 
0.23 
0.20 



40 
50 
50 
50 
50 
40 
10 
40 
40 
10 



52 



20 
36 
56 
58 



22 
22 
24 
28 
28 
22 
22 
22 
36 
36 



22 
20 
24 
20 
20 
20 
22 
22 
36 
28 



19.6 
27.6 
33.2 
35.8 
39.4 
42.4 
44.4 
46.4 
49.4 
53.8 



57 



15 
24 
36 
10 



16.6 
20.6 
26.2 
30.4 
32.0 
36.0 
38.2 
38.0 
46.6 
48.4 



16.6 
20.6 
26.6 
27.0 
29.2 
33.4 
37.8 
38.8 
47.0 
43.8 



Note: a The numbers shown in this column interpolate the results presented in 
Ref. 9 



100 keV, or if the quantity of interest averages the density of states on a typical 
scale of about 100 keV around the Fermi surface. The typical scale which has been 
found, 100 keV, is in fact related to the average inter-distance energy between the 
unbound states which are discretized in the W-S box. In this paper, we clarify the 
effect of the discretization of the unbound states on the ground state energy as 
well as the role of pairing correlations. Consequences for the equation of state are 
presented. 

2. Equation of state in the inner crust 

The inner crust matter is divided in 11 zones shown in Tab. [TJ as in Ref. 4 (notice 
that the denser zone has been removed since it is in the deformed pasta region). Each 
W-S cell is supposed to contain in its center a spherical neutron-rich nucleus sur- 
rounded by unbound neutrons and immersed in a relativistic electron gas uniformly 
distributed inside the cell. For a given baryonic density ps = A/V, neutron number 
N and proton number Z, the total energy of a W-S cell, E to t has contributions com- 
ing from the rest mass of the particles, E m (N, Z) — Z(m p + m e ) + Nm n , from the 
nuclear components (including Coulomb interaction between protons) E n (pB, N, Z), 
from the lattice energy, 



E^pB.N.Z^pp}) = - / d 3 ri d 3 r 2 p e ( ri )- r(p e (r 2 ) -p P {r 2 )) , (1) 



which is induced by the difference between the uniform electron density and proton 
density, and finaly from the kinetic energy T e of the relativistic electron gas. The 
Coulomb exchange energy as well as the screening correction for the electrons gas 
has been neglected. 

The self-consistent Hartree-Fock-Bogoliubov (HFB) approach has been pre- 
sented in detail in Ref. 14 It has been extended to describe the dense part of the 
inner crust at finite temperature in Ref. 15 With minor modifications, it has been 
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Fig. 1. The binding energy e/ps of the W-S cells in zone N ZO ne='l as well as the partial contri- 
butions of the rest mass energy E m /A — m n , the nucleons, the electrons and the lattice. HF and 
HFB results are shown on the left and right panels, respectively. See the text for more details. 



successfully applied to calculate the specific heat of the whole part of the inner 
crust. 16 The effective N-N interaction should provide a reasonable description of 
both the nuclear clusters and the neutron gas. We use the Skyrme force SLy4 17 
which was adjusted to describe properly the mean field properties of neutron-rich 
nuclei and infinite neutron matter. The choice of the pairing force is more problem- 
atic since at present it is not yet clear what is the strength of pairing correlations in 
low density neutron matter. However, in the present paper we have chosen to adjust 
the effective pairing interaction on the results of the pairing gap obtained with the 
bare interaction. We then have a maximum gap in neutron matter of about 3 MeV 
located at Uf — 0.85 fm _1 . 

The minimization is performed in two steps: first, at a given density ps and for a 
given mass number A, we look for the cell which satisfies the beta-stability criterion, 
fip + fj, e = p n . The neutron and proton chemical potentials are provided by the HFB 
calculation while the electron chemical potential is deduced from p e = dE to t/dN e . 
Note that, due to shell effects we may find several cells which satisfy this criterion. 
We then choose the one which has a minimal energy. In the second step, we minimize 
the total energy B to t with respect to the mass number A. It should be noted that, 
in the second step all the cells must be calculated at the same density while the 
minimum energy is searched over the variable A. As the volume is modified by unit 
steps (the step of the mesh is 0.2 fm), the mass number A can take only discrete 
values so that A — psV is always satisfied. The step between two values of A is not 
constant and vary as the square of the WS radius Rws- We have approximatively 
dA ~ p B ^Rws dR ws- 

For N zone —4:, we show in Fig.[T]a typical example of step 2 in the search for the 
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W-S cell which minimizes the energy. On the l.h.s. is represented the total binding 
energy E to t/A of the W-S cells obtained in Hartree-Fock (HF) mean- field model 
(without pairing) . The position of the lowest minimum is indicated as well as two 
other local minima. In this case, there is a difference of 5 keV per nucleon between 
the two first minimum. On the r.h.s., we plot the results of HFB calculations which 
show the effects of the pairing correlations: the first minimum is deeper (about 
34 keV per nucleons). On the lower part of the graph, we show the partial con- 
tributions of the rest mass energy E m /A — m n , the nuclear binding energy E n /A, 
the electron kinetic energy T e /A and lattice energie Ei/A. One sees that the inner 
crust is a typical example of frustrated system: when the nuclear energy is mini- 
mum, the electron energy is maximum and vice-versa. The minimum in energy is a 
compromise between these two opposite energies. The consequence is that it gives 
rise to numerous local minima which are not very far in energy from the absolute 
minimum. 

The equation of state obtained at the HF and HFB approximations is given 
in Table [TJ In total, 19714 cells configurations has been calculated. The pairing 
correlations do not strongly modify the HF results, it may slightly decrease the 
proton number and the WS radius. We also show previous calculations of equations 
of state. We have obtained proton numbers Z different from those of Refs. 4 ' 9 The 
results of Ref. 4 have been obtained in a density matrix expansion framework in- 
terpolating microscopic calculations for homogeneous matter. The difference may 
be due to the fact that the Skyrme interaction SLy4 is based on more recent mi- 
croscopic parametrizations. It is difficult to compare with the results of Ref. 9 since 
those calculations have not been done at the same densities and they have essen- 
tially been carried out in the dense part of the crust. We then show an interpolation 
of their results. The comparison shows a similar trend in the dense cells. 

3. Discussion of the shell effects 

In ordinary nuclei, the bound states play an essential role in determining the ground 
state properties and the role of the continuum states is usually negligeable. The cou- 
pling between bound states and continuum states can be more important for nuclei 
near the drip lines, as the Fermi energy is increasing. However, even in these most 
extreme situations continuum states are only virtually populated and the occupa- 
tion numbers associated to continuum states are small compared to those of the 
bound states. In the inner crust of neutron stars, the situation is completely dif- 
ferent: the Fermi energy is positive and most of the mass is due to the unbound 
neutrons while bound states account only for a small fraction. Thus, the continuum 
states contribute in an essential way to physical properties of the crust like its den- 
sity, total energy or specific heat. The presence of bound states, or nuclear clusters, 
is not negligible but it induces only corrections to the ground-state properties. 16 It 
is then very important to have an accurate model of the unbound states. 

It should also be noted that, as the continuum states are embedded in a Coulomb 
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lattice, their asymptotic behavior is different from that occurring in nuclei. The 
continuum states are indeed periodically distorted within a distance equal to the 
lattice step (20-100 fm.). The proper theory for the description of the continuum 
states is the band theory which requires to solve the Schrodinger equation, 12,13 
(ho + /ik)w Q ,k(r) = £a,kW Qi k(r), where ho is the usual single-particle Hamiltonian of 
mean-field models. The second term is induced by the Floquct-Bloch boundary 
conditions for the single-particle wave function ip a y. <^ a ,k( r ) = Ua,k(r)e lk r . The 
spherical W-S approximation is obtained by setting k = and replacing the W- 
S elementary polyhedron by a sphere. Then, the remaining index a is discrete, 
like in usual mean-field models, and runs over the eigenstates of the Hamiltonian 
ho- The main difference between mean- field models and band theory is indeed due 
to the term h^. In band theory, this term introduces a new index k which is a 
continuous variable ranging between and n/a fm -1 where a is the lattice step. As 
a consequence, the density of unbound states is continuous (notice that it may also 
have structures and gaps). It is then clear that the discrete distribution of unbound 
states present in usual mean-field models is a consequence of the dropping of /ik in 
the W-S approximation. 

We can characterize the distribution of unbound states by the average distance 
between the unbound states energies, which can be related to the radius of the W-S 
cell as Ae ~ h/2mR^ vs : the smaller the W-S radius, the larger the average distance 
Ae. One could notice that the W-S radius is a function of the density and it is 
fixed by the equation of state, as shown in Table [T] At the W-S approximation, the 
distribution of unbound states is then also fixed by the density. The non-continuous 
distribution of unbound states could be interpreted as a spurious effect, induced by 
the W-S approximation, since in band theory calculations, the density of unbound 
states varies in a continuous way. 13 One can then distinguish between the physical 
shell effects due to the bound states from the spurious shell effects due to the 
unbound states and the W-S approximation. 

4. Spurious shell effects 

In order to estimate the spurious shell effects due to the unbound neutrons, we 
simulate an homogeneous gas of neutron matter in W-S cell by removing all the 
proton states. In homogeneous matter, the energy density should be independent 
of the volume, then of the W-S radius, and it should depend only on the density. 
In our simulation of homogeneous matter, we then vary the W-S radius from 10 
to 50 fm, at constant neutron density. We call BwS— homXPunb., Rws) the binding 
energy per particle obtained for a given W-S radius Rws and for a given density 
of unbound neutrons p un b. ■ We show the binding energy versus the W-S radius in 
Fig.[2^a) for several densities corresponding to N zone —6, 7 and 8. It is clearly seen 
that the binding energy is not independent of the W-S radius, except for large W- 
S radii where the binding energy converge to the energy of homogeneous matter, 
BhomXPunb.)- The binding energy of homogeneous matter BhomXPunb.) is calculated 
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Fig. 2. We represent several quantities versus the W-S radius: (a) the binding energy 
BwS-hom.{punb.,Rws) corresponding to 7V zone =6, 7 and 8 (b) the difference B hom \p mb ) - 
B\v s—hom.(Punb.i Rw s) (dots) and the function f(p U nb.,R\vs) which fit the smooth compo- 
nent (solid lines), and (c) the difference between the dots presented in (b) and the function 

fiPunb., Rws)- 

from the average value of the Skyrme two-body interaction on a plane wave basis. 
The difference B hom {p unh ) - B ws _ hom , {p m b. , Rws) is represented on Fig. [2b) 
(dots) . As expected, this quantity is converging to zero for large W-S radii, but for 
small radii, around 20 fm for instance, the difference can be as large as 300 keV per 
nucleons. 

From the dots represented on Fig. E^b), we obtain an universal func- 
tion which fits the smooth radial dependence of the difference Bhom.{Punb.) — 
Bws-hom.{Punb., Rws)- As this effect is due to the discretization of the unbound 
states, it is proportional to the average distance Ae between the box states. This 
difference should then decrease with the W-S radius as R^g- It should also be pro- 
portional to the number of unbound neutrons, while the density dependence is not 
known. We then represent it by a power law p" nb where the power a is adjusted to 
reproduce the dots presented in Fig.[2]Jb). In this figure we show the fitted function, 

f(p mb ,Rws) = sg.os^./po) 1425 ^ • ( 2 ) 

It interpolates the smooth part of the shell effects on the binding energy. There 
is, however, a residual difference which cannot be fitted. Finally, to estimate the 
residual difference between Bh om .(Punb.) — Bws-hom.(Punb., Rws) and the adjusted 
function f(p U nb.,Rws), shown on the r.h.s of Fig. [2jc). The average residual fluc- 
tuations are now of the order of 50 keV. 
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5. Modified HF binding energy in W-S cells (HF-WS) 

In Sec. |3J it has been shown that the fluctuations in the binding energy induced 
by the spurious shell effect can be decomposed into a smooth and a residual term 
and we have obtained a universal function for the smooth term. We then propose a 
systematic correction to the HF binding energy in the W-S cell: 

BhF-Ws{PB, Punb., Rws) = By/ S — irihom. (pb,Rws) 

+ ( Bhom. (Punb. ) — BwS-hom. {Punb. , Rws) ) 

V / smooth 

= B ws — irihom. (PB, Rws) + f (punb., R\VS) ■ (3) 

In the calculation of the total energy, the nuclear binding energy E n /A should then 
be replaced by Bhf-ws{pb, Punb., Rws)- The minimum energy can be unambigu- 
ously identified if the difference in energy between the first and the second minimum 
is larger that the residual fluctuation, about 50 keV. In the HF calculation, the dif- 
ference in energy between different cell configurations are of the order of few keV 
to few tens of keV. The residual fluctuation of the proposed correction is then too 
large. However, as seen above, the pairing correlations help in producing a deeper 
first minimum. It may also help in reducing the size of the residual fluctuations. In 
a future investigation, we plan to obtain an improved binding energy for W-S cells 
including pairing correlations. It must be remarked that the effects of the coupling 
between the protons and the unbound neutrons have been neglected in Bhf-ws- 
This is a good approximation since the protons are deeply bound. 

6. Conclusions 

In the search for the W-S cells which minimize the energy, it is important to look not 
only at the lowest minimum energy, but also at the other minima. The differences 
between the first and other minima show how stable is the lowest minimum with 
respect, for instance, to thermal fluctuations or spurious shell effects. In HF calcula- 
tions, the differences in total energy between the first minimum and the others are 
about few keV per nucleon (few MeV in total energy). It is shown that the spuri- 
ous shell effects, induced by the W-S approximation, introduce a fluctuation in the 
calculation of the total energy up to about 300 keV per nucleon. Interpolating the 
smooth part of the spurious shell effects, we have then proposed a simple method to 
reduce the fluctuations down to about 50 keV per nucleon (HF-WS). The residual 
fluctuations are still too large, and the removal of the smooth part is not enough to 
improve the results obtained within the W-S approximation. We have also shown 
that the pairing correlations help in producing a deeper first minimum. We have 
then compared the equation of state obtained from HF and HFB binding energies. 
Several issues should still be addressed in future studies: spurious shell correction 
with pairing correlations as well as effects of the temperature. These two effects 
should give lower residual fluctuations than the one we obtained at the T = HF 
approximation. 
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